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Abstract 

We present the applications of methods from nonlinear 
local harmonic analysis in variational framework to cal- 
culations of nonlinear motions in polynomial/rational ap- 
proximations (up to any order) of arbitrary n-pole fields. 
Our approach is based on the methods provided possibil- 
ity to work with dynamical beam/particle localization in 
phase space, which gives representions via exact nonlinear 
high-localized eigenmodes expansions and allows to con- 
trol contribution to motion from each scale of underlying 
multiscale structure. 

1 INTRODUCTION 

In this paper we consider the applications of a new nu- 
merical-analytical technique based on the methods of local 
nonlinear harmonic analysis or wavelet analysis to the cal- 
culations of orbital motions in arbitrary n-pole fields. Our 
main examples are motion in transverse plane for a sin- 
gle particle in a circular magnetic lattice in case when we 
take into account multipolar expansion up to an arbitrary 
finite number and particle motion in storage rings. We re- 
duce initial dynamical problem to the finite number of stan- 
dard algebraical problems and represent all dynamical vari- 
ables as expansion in the bases of maximally localized in 
phase space functions. Our approach in this paper is based 
on the generalization of variational-wavelet approach from 
[1]-[12]. Starting in part 2 from Hamiltonians of orbital 
motion in magnetic lattice with additional kicks terms and 
rational approximation of classical motion in storage rings, 
we consider in part 3 variational-biorthogonal formulation 
for dynamical system with rational nonlinearities and con- 
struct explicit representation for all dynamical variables as 
expansions in nonlinear high-localized eigenmodes. 

2 MOTION IN THE MULTIPOLAR 
FIELDS 

The magnetic vector potential of a magnet with 2n poles 
in Cartesian coordinates is 



A = ^2K n f n (x,y), 



(1) 



where /„ is a homogeneous function of x and y of order n. 
The cases n = 2 to n = 5 correspond to low-order mul- 
tipoles: quadrupole, sextupole, octupole, decapole. The 
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corresponding Hamiltonian is ([13] for designation): 
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Then we may take into account an arbitrary but finite num- 
ber of terms in expansion of RHS of Hamiltonian (||) and 
from our point of view the corresponding Hamiltonian 
equations of motions are not more than nonlinear ordinary 
differential equations with polynomial nonlinearities and 
variable coefficients. Also we may add the terms corre- 
sponding to kick type contributions of rf-cavity: 
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or localized cavity V(s) — Vo-S p (s— sq) with S p (s — sq) = 

X)n=-oj — ( s ° + n ' L)) at position sr> We consider, 
as the second example, the particle motion in storage rings 
in standard approach, which is based on consideration in 
[13]. Starting from Hamiltonian, which described clas- 
sical dynamics in storage rings and using Serret-Frenet 
parametrization, we have after standard manipulations with 
truncation of power series expansion of square root the fol- 
lowing approximated (up to octupoles) Hamiltonian for or- 
bital motion in machine coordinates: 
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Then we use series expansion of function f(p a ) from [13]: 
f(p a ) - /(0) + /'(OK + /" (0)^/2 + . . . = p a - 
p 2 r /(2j 2 ) + . . . and the corresponding expansion of RHS 
of equations corresponding to (4). In the following we take 
into account only arbitrary polynomial/rational (in terms of 
dynamical variables) expressions. 



3 VARIATIONAL APPROACH IN 
BIORTHOGONAL WAVELET BASES 

The first main part of our consideration is some vari- 
ational approach to these problems, which reduce initial 
problem to the problem of solution of functional equations 
at the first stage and some algebraical problems at the sec- 
ond stage. Multiresolution expansion is the second main 
part of our construction. As a result the solution is param- 
eterized by solutions of two reduced algebraical problems, 
one is nonlinear and others are linear problems obtained 
from wavelet constructions and represented as expansion 
in a compactly supported wavelet basis. Because inte- 
grand of variational functionals is represented by bilinear 
form (scalar product) it seems more reasonable to consider 
wavelet constructions which take into account all advan- 
tages of this structure. Let (M, us) be a symplectic man- 
ifold, H is Hamiltonian, Xh is unique Hamiltonian vec- 
tor field defined by uj(Xh(x), v) = —dH(x)(v), v E 
T X M, x E M, where u> is the symplectic structure. T- 
periodic solution x(t) of the Hamiltonian equations x — 
Xh{x) on M is a solution, satisfying the boundary condi- 
tions x{T) = x(0),T > 0. Let us define a function 



i < -Jx,x >dt- I U.(.r{ (5s 



The critical points of $ are the periodic solutions of x — 
Xh{x). Computing the derivative at x E ft in the direction 
of y E CI, we find 

&(x)(y) = ^rHx + ey)\ e =o (6) 
de 

= < -Jx - xjH{x), y > dt 
Jo 

Consequently, $'(x)(y) = for all y E fi iff the loop x 
satisfies the equation 



Jx(t) - S7H(x(t)) = 



(7) 



Now we introduce wavelets in our underlying bilinear 
structure. Let us consider action of operator S on x 

S(cj(J),H, x, d/dt, V, t)x = ~Jx{t) - xjH{x{t)) (8) 

which is polynomial/rational in x, and have arbitrary de- 
pendence on t. Then (6) is equivalent to 
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We start with two hierarchical sequences of approximations 
spaces [14]: 

...V- 2 cV- 1 cVoCV 1 cV 2 ..., (10) 
...V-2 CV-i C Vq C Vi CV 2 ..., 
and corresponding biorthogonal expansions: 
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Let Wq be complement to Vq in V\, but not necessarily 
orthogonal complement. Orthogonality conditions have the 
following form: W 1 Vq, Wq 1 Vq, Vj _L W 3 , Vj ± 
Wj . Translates of ip span Wq, translates of tp span Wq. 
Biorthogonality conditions are 



< ipjk,i>3>k> >= / ip 3 k{x)ipyk'(x)Ax = S k k'Sjj', 
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(12) 

where ^jk(x) — 2^ 2 %p{2^ x — k). Functions tp(x), (p{x—k) 
form dual pair: < <p(x — k),tp(x — £) >= Ski, < <p{x — 
k),ip(x — £) >= 0. Functions ip, (p generate a multiresolu- 
tion analysis. ip(x — k), ip(x — k) are synthesis functions, 
(p(x — £), ip(x — £) are analysis functions. Synthesis func- 
tions are biorthogonal to analysis functions. Scaling spaces 
are orthogonal to dual wavelet spaces. Two multiresolu- 



tions are intertwining Vj+Wj = Vj+i , Vj +Wj 
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These are direct sums but not orthogonal sums. So, our 
representation (11) for solution on the level of resolution 



Vj has now the form 



(13) 



where synthesis wavelets are used to synthesize the func- 
tion. But djk come from inner products with analysis 
wavelets. Biorthogonality yields 



x j (t)i>j m (t)dt. 



(14) 



So, we may use this more useful construction in our vari- 
ational approach [1]-[12]. We have modification only on 
the level of computing coefficients of reduced nonlinear al- 
gebraical system of equations. This biorthogonal construc- 
tion is more flexible and stable under the action of large 
class of operators while orthogonal (one scale for multires- 
olution) is fragile, all computations are much more simpler 
and we accelerate the rate of convergence. In all types of 
Hamiltonian calculation, which are based on some bilin- 
ear structures (symplectic or Poissonian structures, bilinear 
form of integrand in variational integral) this framework 
leads to greater success. In numerical modelling we may 
consider very useful wavelet packets. As a result we have 
from (9) the following reduced system of algebraical equa- 
tions (RSAE) on the set of unknown coefficients of ex- 
pansions (11): 



L(Sij,a,ai,f3j) = 
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where operator L is algebraization of initial problem (8). 
/ = ...,i q+2 ), J = (ii, jp+i) are multiindexes, 
by which are labelled ai and the other coefficients of 
RSAE (|15|): 



■ jp+i } J 



n vs* 

i<jfc<p+i 



(16) 



0:1 



(f il ...ip is ...ip iq+2 , 



where p (q) is the degree of nominator (denominator) part 
of operator S (8), ii = (1, q + 2), ip'j H = d<pi B /dt. Now, 
when we solve RSAE ( |l5| ) and determine unknown coeffi- 
cients from formal expansion (1 1) we therefore obtain the 
solution of our initial problem. It should be noted that if 
we consider only truncated expansion with N terms then 
we have from ( |i"5| ) the system of N x n (n is dimension 
of x (7)) algebraical equations and the degree of this alge- 
braical system coincides with degree of initial differential 
system. The problem of computations of coefficients ai, 
Pj ( |l6| ) of reduced algebraical system may be explicitly 
solved in wavelet approach. The bases functions ^(t) (11) 
are obtained via multiresolution expansions (10) and repre- 
sented by compactly supported wavelets. Because affine 
group of translations and dilations is inside the approach, 
this method resembles the action of a microscope. We have 
contribution to final result from each scale of resolution 
from the whole infinite scale of spaces (10). The solution 
has the following form 

x(t) = xfj ow (t) + x i M)' u i ~ 23 

j>N 

which corresponds to the full multiresolution expansion in 





Figure 1: Multiscale representations for x — p x sections. 

all time scales. Formula (|l7|) gives us expansion into a slow 
part zff w and fast oscillating parts for arbitrary N. So, we 
may move from coarse scales of resolution to the finest one 
for obtaining more detailed information about our dynam- 
ical process. The first term in the RHS of equation (19) 
corresponds on the global level of function space decom- 
position to resolution space and the second one to detail 
space. In this way we give contribution to our full solution 
from each scale of resolution or each time scale. On Fig. 1 



we present sections x — p x corresponding to model (5) in 
different parameter regions. 
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